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qq , We propose a simple parametrization for the deep-inelastic diffractive cross section. It con- 

tains the contribution of qq production to both the longitudinal and the transverse diffractive 
structure functions, and of the production of qqg final states from transverse photons. We start 
from the hard region and perform a suitable extrapolation into the soft region. We test our 
model on the 1994 ZEUS and HI data, and confront it with the HI conjecture of a singular 
gluon distribution. 
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1 Introduction 



Diffractive events are characterized, in general, by the presence of large rapidity gaps in the hadronic 
final state that are not exponentially suppressed. These are conventionally ascribed to Pomeron 
exchange. Diffractive processes in deep-inelastic scattering (DIS) are of particular interest, because 
the hard photon in the initial state gives rise to the hope that, at least in part, the scattering 
amplitude can be calculated in perturbative QCD (pQCD). With the increasing amount of data 
on diffractive DIS [jl], ||] , we have reached a level of accuracy that provides deeper insight into the 
nature of the Pomeron and its coupling to partons. Rapidity-gap events make up a sizable fraction 
of all DIS events, and can only be due to the exchange of some colourless object in the cross channel. 
The simplest realization of the Pomeron in pQCD is provided by two gluons of opposite color (3|. 
More detailed models based upon the two-gluon picture have been formulated both for perturbative 
gluons [|j and for nonperturbative (massive) gluons ||, §, 0, |8|: for an alternative approach, see ||. 
In this paper we take the point of view that perturbative QCD provides a reasonable starting 
point, as in the more detailed model described in [||, [Hi. In the present paper, we develop a simple 
parametrization based on this model, which can easily be compared with experimental data. 

The simplest description of diffractive DIS starts from the process of qq production. This 
process is conveniently described in terms of light-cone wave functions [§,11, 0] in the proton rest 



frame. The light-cone wave function of the photon contains the information about the dissociation 
of the fast-traveling photon into partons, long before the interaction with the proton occurs. At the 
beginning of the scattering process, the photon splits into a quark-antiquark pair. At sufficiently 
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large photon virtuality Q 2 , the quark-antiquark pair radiates additional gluons before it reaches 
the proton at rest. At the time of the interaction, the partonic system is spread over a transverse 
area which is comparable with the size of the hadron. One expects therefore that the exchanged 
Pomeron should be close to the usual soft hadronic Pomeron. However, inside the final state of 
the partonic system, we expect that there are also "hard" configurations, for which the exchanged 
Pomeron should behave quite differently. These are final states for which the partonic system is 
confined to small transverse distances. Examples are longitudinally-polarized vector particles [12] 



and high-pr jets JT^, 14]. In the inclusive measurement of diffractive final states, one sums over 
both these small-distance and large-distance configurations. So far there is no theoretical framework 
which allows one to predict the relative magnitudes of the "soft" and the "hard" components of 
the diffractive cross section, which must be determined by experiment Since the cross section for 
the "soft" component is expected to rise weakly with energy for any fixed mass of the diffractive 
system, whereas the "hard" part should rise faster, the energy dependence of the diffractive cross 
section may help to determine the relative sizes of the two components. 

In attempting to formulate a model that interpolates between these two components, one finds 
that perturbative models based upon two-gluon exchange, which are valid a priori only for small- 
size final states, allow a smooth extrapolation into the soft region. Thanks to gauge invariance and 
colour cancellation, one does not encounter infrared singularities in models for qq production or 
qqg production, i.e., there is no need for an artificial cutoff. Moreover, the wave- function formalism 
can be extended to include multi-gluon exchange, which is useful for going beyond the lowest-order 
two-gluon exchange. 

Since the first observation of diffractive DIS at HERA, several attempts have been made to 
compare the data with QCD-based models. In particular, the concept of the Pomeron structure 
function and its DGLAP Q 2 evolution has been applied 16, 17, [l8|, gj. In these analyses, the 
diffractive cross section has been assumed to consist only of leading twist, and the contribution of 
the longitudinal photon has been disregarded. On the other hand, there is strong evidence that, for 
small masses of the diffractive system, the longitudinal cross section for qq production - although 
formally of higher twist - is not small compared to the transverse cross section. We therefore feel 
that a more complete analysis of the HERA data should include the longitudinal cross section, 
and particular emphasis should be given to the region of small diffractive masses. As a minimal 
model, one might consider just the production of quark-antiquark pairs, which should dominate 
the small-mass region. For somewhat larger masses, the production of an extra gluon has to be 
taken into account. In a future step, one would also have to address the Q 2 evolution of both the 
transverse and the longitudinal cross sections. 

In this paper we propose and test a simple parametrization of the diffractive cross section that 
is motivated by the above considerations. Stimulated by perturbative QCD, we make an ansatz 
that consists of four terms, which model both the transverse and the longitudinal cross sections for 
qq and qqg production. By treating the overall strength of these terms and the exponent of the 
energy dependence as a free parameter, we let the data decide which fractions of the diffractive 
cross section belong to the "soft" and "hard" parts. After a brief description of the model, we 
compare this Ansatz with the ZEUS and HI 1994 data, and finally draw a few conclusions. 



1 It was pointed out in [[L5| that the pseudo-rapidity cuts imposed in certain early analyses selected a "hard" 
component in diffractive DIS, which must also be present, at some level, even in analyses without this a priori 
selection. 
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2 A Parametrization for Diffractive DIS 
2.1 Theoretical Motivation 

The main variables used for the description of diffractive DIS are the total hadronic energy W 
of the 7*-proton system and the diffractively-produced mass M. In the analysis of the diffractive 
structure function, it is convenient to use also the variables (3 and xjp. In terms of W and M, one 
has (3 = Q 2 / (M 2 + Q 2 ) and xp = (M 2 + Q 2 ) / (W 2 + Q 2 ), where we have neglected the proton mass 
and the momentum transfer t. To connect these variables with the Bjorken scaling variable xb, we 
recall that xb = Q 2 /(W 2 + Q 2 ), which immediately leads to xb = (3xjp. 

Before describing our model in somewhat more detail, we first make a few general remarks. 
First, we expect the cross section to be dominated by very small t values. After integration over 
final-state kinematic variables, the t dependence and the strength of the coupling of the Pomeron to 
the proton will be combined in the overall normalization. Next, the j3 spectrum and the (5 2 -scaling 
behavior follow from evolution of the final-state partons, and can most easily be derived from the 
light-cone wave functions of the incoming photon. This part of our model therefore decouples from 
the dynamics inside the Pomeron. On the other hand, the energy dependence of the Pomeron, 
i.e., the xjp-distribution, can be calculated at most only partially within perturbative QCD. We 
therefore leave it free parameter. 

We have already indicated that, in the proton rest frame, the light-cone wave function formalism 
provides a nice intuitive description of diffractive DIS. At leading order, when the photon dissociates 
into a quark-antiquark pair, we have a single color dipole with a certain momentum distribution 
given by the corresponding wave function. At higher order, more partons are generated and the 
initial state can be rather complex. At leading-twist level, however, the basic structure is again 
a single color dipole: all partons (gluons or quarks) but one are located within a small area in 
impact-parameter space, i.e., at short relative distances, whereas the remaining single parton is 
well separated. The localized parton subsystem carries color conjugate to that of the single parton, 
so that one has again a color-dipole configuration. The short-distance evolution within the parton 
subsystem factorizes, so one only needs to introduce a wave function for the momentum distribution 
of the single parton. With this simplification, we end up with two basic structures: a quark- 
antiquark dipole and a gluon-gluon dipole, the latter appearing only at higher order. 

In the case of the elementary quark-antiquark final state, the wave function depends on the 
helicities of the photon and of the (anti)quark. We define left- and right-handed transverse photons 
by projecting on the polarization vectors (l,i) and (1,-i) (7 = ±1), respectively, and the longitudinal 
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polarization vector is proportional to the proton momentum p (7 = 0). For massless quarks, the 
spin is orientated along the direction of motion or opposite to it (h = ±1). As variables for the 
wave function, we use the Sudakov parameters k = aq' + P^P + kt with q' = q + %bP- I n the 
proton rest frame with a fast-moving photon, the parameter a is of the order unity and denotes 
the momentum fraction of the photon momentum carried by the quark (Figs, la, 2), whereas /9& 
is small and may be neglected. Using a complex notation for kt- k = k x + ik v , k* = k x 



finds for the transverse photon (see also [11]): 



iky, one 



V2 (1 - a) k t 
\k t \ 2 + a(l - a)Q 2 

y/2 a kt 
\k t \ 2 + a(l - a)Q 2 

V2akt 
\k t \ 2 + a(l - a)Q 2 



for 7 = +1 and h = +1 
for 7 = +1 and h = — 1 
for 7 = —1 and h = +1 



(1) 



V2(l 



a 



kt 



\k t \ 2 + a(l - a)Q 2 



for 7 



-1 and h 



-1 



Similarly, for the longitudinal photon one finds: 



^l(a,k t 



a(l — a) Q 
\k t \ 2 + a(l - a)Q 2 



for 7=0 and h = ±1 



(2) 



Here Q 2 , 7 and h denote the virtuality of the photon and the helicities of the photon and the 
quarks, respectively, and the wave function includes the propagator of the off-shell quark carrying 
the momentum k. In the following, we shall make use of the small-fej behavior of (1) and (2). 

To produce a gluon dipole (Fig. lb), we have again to start from the approximation described 
previously, because a direct coupling of photons to gluons is lacking. In the leading-log (Q 2 ) ap- 
proximation, transverse momenta are strongly ordered, which translates into the inverse ordering 
of distances in impact-parameter space. The quark-antiquark pair with a large transverse mo- 
mentum is localized in impact-parameter space, and forms an effective "gluon" state conjugate in 
color to the emitted gluon, which has a smaller transverse momentum and is separated by a large 
distance from the quark-antiquark pair. In this approach only the transverse photon polarization 
is of importance, and it determines the hard part of the process, i.e., the effective gluon dipole is 
independent of whether the photon is right- or left-handed. The wave function of the gluon dipole 
has the following tensor structure (the indices fi,v = I, 2, since only the transverse components are 
involved) : 

1 h 2 n^ v - 1 k^k v 

W(a h) = 1 9t 1 1 (3) 

{ ' t} y/a(l - a)Q* k 2 + a(l-a)Q 2 ' W 

where a and kt refer to the diffractively-produced gluon. We have written explicitly the (1 — a) 
term, even though the approximation applied here requires a to be much less than 1, so that 
(1 — a) ~ 1. We have introduced this term in order to make manifest the analogy to the previous 
expression for the quark dipole. Our main interest is, again, the behavior of (3) near kt = 0. 

So far we have discussed the photon wave functions which describe the qq or qqg state after the 
splitting of the photon into the quark-antiquark pair. In order to obtain the diffractive scattering 
amplitude, we have to include the interaction with the proton target. Beginning with the qq final 
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k+xpp 



a b 
Figure 2: Interaction with the proton, modelled in the two-gluon exchange approximation. 



state, this interaction is represented in Fig. 2. It is essential that the two gluons couple in all 
possible ways to the two quarks. In the proton rest frame, the lower quark with momentum k emits 
a gluon which, after the interaction with the proton, is reabsorbed by one of the two quarks. A 



more detailed discussion of this two-gluon-exchange model has been given elsewhere [10]. Here we 
only briefly describe a few main features that we need in order to motivate our parametrization. 
The coupling of two i-channel gluons with zero net color and transverse momenta It and —It to the 
color dipole, i.e., to any wave function of type (1), (2), or (3), can be obtained simply by taking 
differences of the wave function: 



1 t dkldkl 1 



2*(a,fct) 



for l t 







oo 



(4) 



and then convoluting with a suitable ansatz for the Zf dependence of the Pomeron form factor of 
the proton. 

It is an important feature of our two-gluon model that our simple wave functions (1), (2), 
(4), together with a suitable ansatz for the Pomeron amplitude, provides an interpolation between 
the hard and the soft region. In particular, one finds that the transverse polarization belongs to 
leading twist and is dominated by the aligned-jet configuration, whereas the longitudinal polar- 
ization of the photon leads to a higher-twist contribution. In order to see this behavior in the 
wave-function formalism, we note the relation between the diffractive structure functions F D and 
our wave functions 



F D ( XjP ,f3,Q 2 ) ~ Jdt J 



l 2 t 



D^(a,k t ) </>(!?,k%;x P ) 



(5) 



where D^> is taken from (4), and (f>(l 2 , k^xjp) stands for the Pomeron amplitude. Here fcg denotes 
some hadronic scale which separates the regions of soft and hard QCD: it should not be confused 
with the QCD factorization scale. The variables a, (3, and M 2 are related through 



a(l - a)M 2 



k 



t • 



(6) 
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and (3 = Q 2 /{Q 2 + M 2 ). A simple choice for the It dependence is [O] 

(7) 

where K^Ao) k 1 as I 2 > fcg, and v{l 2 /kl) -> as / 2 -> 0. 

We start in that part of the kinematic region where the parton model of Fig. 2 is most reliable, 
i.e., the region where the virtuality of the quark with momentum k is large: k 2 + a(l — a)Q 2 = 




k 2 /(l — (3) > fcg [13]. In this final-state configuration, the g<j pair has a small transverse size, 
and the two-gluon Pomeron interacts with the whole system. Consequently, both contributions of 
Fig. 2 are important, and the simple picture of a "Pomeron structure function", which would be 
suggested if only Fig. 2a were taken into account, does not apply 0. In this region, the Pomeron 
amplitude (j)(l 2 , k 2 ,; xjp) coincides, to a good approximation, with the unintegrated gluon structure 
function of the proton: 

r*?/(l-0 



/ 



dlt<P(lt,k^ xjp) = xjpgixjp, kt/(l - (3)) (8) 



(for the extrapolation into the region of smaller kt and for a discussion of the dependence on xp, 
see below). Going into more detail, let us look into the dependence upon k 2 and a at fixed Q 2 . In 
the It integral, it suffices to note that the dominant region is kfi < I 2 < k 2 + a(l — a)Q 2 : in our 
example (7), <f> ~ l// 2 , and the dominant contribution, in fact comes from this kinematic domain. 
In this region, we approximate D^f in (4) by the limit It — > 0, and obtain: 

k 2 +a{1 _ a)Q 2 2l a(l-a)\k t \Q 2 , )2 M n ^ 2n 

Inserting this into (5) and making use of relation (6), we find that the integral over k 2 (at fixed 
j3) is dominated by the lower limit k^. In terms of the variable a, this lower limit corresponds 
to a ~ k\jQ 2 or 1 — a ~ k^/Q 2 . The final result for (5) is constant in Q 2 , i.e., it is of leading 
twist. The end points of the a integral correspond to the aligned configuration: in the center-of- 
mass system of the quark-antiquark pair (Fig. 3), a is related to the scattering angle 6 through 
2a = 1 — cos 9, and the dominant regions are 6 = 0, n. In other words, starting in the hard region 
of large transverse momenta, we find that the main contribution comes from the lower end of the 
kt integral, i.e. we find ourselves pushed into the soft region where our perturbative ansatz for the 
Pomeron becomes invalid. 

A similar argument applied to the longitudinal case shows that the dominance of small k 2 (or 
values of a close to zero or one) is less pronounced: instead of (9), we now have 

k? +a (i- a) Qi d 2 k ^ 2 _ a(l-a)VQ^(a(l-a)Q 2 -k 2 ) ______ u2 



j2 </>(l t ,k ;x F )D^ ~ (fc f 2 + Q !(l-a)Q2)3 x F g(x F ,k t + a{l - a)Q' 



„2 



xjpg{xjp,k 2 /{l- f3)). (10) 



kfQ 

Inserting this into (5), we see that for small k 2 the integral diverges only logarithmically. Also, in 
contrast to the transverse case, the result is of order 1/Q 2 and hence belongs to nonleading twist. 

2 This point has previously been emphasized in [^5). 
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Figure 3: Two-jet production in the 7*-iP center-of-mass system 



On the other hand, the integration over k 2 now yields an additional logarithm in Q 2 / (4/3/cq), which 
is absent in the transverse leading-twist case, and slightly compensates for the Q 2 suppression. 

Next we return to the transverse case and take a closer look at the soft region where k 2 + a(l — 
a)Q 2 < k 2 ,, i.e., k 2 < (1-/3)^ and a < j3k 2 /Q 2 or 1 — a < fik^jQ 2 . Now the quark with momentum 
k in Fig. 2a, before it interacts with the two-gluon Pomeron, is nearly on shell, and one expects 
the picture of the Pomeron structure function to become valid, i.e., the lower parton in Fig. 2a 
can be considered more as a "valence" constituent of the Pomeron, and the contribution of Fig. 2b 
should be less important. This transition can be made explicit by changing the I 2 dependence of the 
Pomeron amplitude, which now can no longer be identified with the unintegrated gluon structure 
function, in such a way that is gives more weight to the region I 2 < k 2 ,: the simple example in 
eq.(6) leads to <fi ~ l/^o = constant. The correct behavior is obtained simply by taking in (4) I 2 
much larger than k 2 , i.e., It — ► oo: in this region the first term 2 ^>{a,kt) dominates, and we end 
up with 

dHt 0(l 2 )D* ~ ;2 „ 2 lh } — -. \J e ° ^ 



1} k 2 (k 2 t + a{l - a)Q 2 )' \k 2 + a(l-a)Q 2 




k§ fot 



(11) 



which corresponds to the planar diagram Fig. 2a (on the rhs of (11), we have disregarded the 
dependence upon xp). Returning to (5), the integral in k 2 can easily be performed and leads to a 
finite leading-twist result. By similar arguments the longitudinal case is found to be of the order 
1/Q 2 . 

In summary, with a suitable Ansatz for the It dependence of the two-gluon Pomeron or, even 
more simply, with a simple prescription for the It integral, it is possible to interpolate between the 
hard region where the parton model applies and the soft region where the aligned jet configuration 



dominates. For the latter case we seem to have arrived at the same conclusions as Bjorken [21], 
although with a somewhat different line of arguments. But, as we shall argue further below, there 
is a new element that we have to take into account, namely the observation of the strong rise of the 
gluon structure function at small x, which gives a weight to the region of large transverse momenta 
which is larger than was anticipated before the advent of the HERA data in, for example, ||21| . 

Before we turn to the xjp dependence of the cross section, let us mention that the momentum 
dependence of the qqg final state is quite analogous to that of the qq system. We do not show 
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the analogue of Fig. 2: it is again essential that the two gluons couple to the diffractive system in 
all possible ways. As before, one has to start in the region where the transverse momenta of all 
three partons are large. When trying to integrate over the transverse momentum of the gluon, one 
finds dominance by the low-momentum region: in this region, all nonplanar couplings of the two 
gluon lines of the Pomeron to the diffractive state become less important, and we are left with 
the leading-twist 'Pomeron structure function' picture, where the Pomeron interacts only with the 
gluon and not the quarks. 

Next we adress the xjp dependence of the cross section. So far we have drawn a rather simple 
picture of the diffractive final state: both for the qq and the qqg final state, we have argued that 
in the preferred configuration at least one of the final-state partons has a rather soft transverse 
momentum, and it is this parton which couples to the Pomeron. If the virtuality of this parton is 
characterized by a typical hadronic scale ~ Aqcd, this seems to imply that the energy dependence 
of the diffractive cross section should be the same as in hadron-hadron scattering, i.e., the diffractive 
structure function F® grows as (l/xjp) njp with njp = 2a/p(0) — 1 « 1.12. However, because of the 
observed rise of the gluon structure function at small x, the situation is more complicated. Let 
us return to the above discussion of the qq final state. The perturbative region is that of large 
transverse momenta of the final-state partons. For this part of the phase space, we expect the 
Pomeron to be described by the perturbative two-gluon model, i.e., the x dependence of the cross 
section will be given by the square of the gluon structure function of the proton (ij 



J jL ml „ [XFg{XF y t , (1 -( i))? {12) 

This should lead to a rise ~ (l/xip) nhard where nhard = lot-hard, — 1 grows with the transverse 
momentum k 2 of the partons, and typically lies above the value 1.4. For the kinematic region where 
the quark transverse momenta are small and our perturbative two-gluon Pomeron has to be replaced 
by some model for the nonperturbative Pomeron, we expect a smaller exponent n: the conventional 
soft Pomeron would suggest that n = njp = 2ajp(0) — 1 ~ 1.12. Since in the diffractive cross section 
we integrate over both the perturbative and nonperturbative parts of the phase space, there will 
be competition between the two regions. At first sight, the large-momentum region seemed to be 
rather subdominant. However, the large gluon structure function provides an enhancement of this 
region, and in this way weakens the dominance of the soft nonperturbative region. As a result, the 
effective value of the exponent n, n e ff, is expected to lie somewhere between the soft and the hard 
values, and the effective scale at which the kt integral peaks should be somewhat higher than the 
soft Pomeron scale. Theoretical studies [19| indicate that n e // only weakly depends upon Q 2 , but 
they do not allow us to predict the numerical value of n e ff or the momentum scale. 



2.2 The Parametrization 

After this brief theoretical review, we are ready to describe and motivate our parametrization. It 
will be given in terms of the diffractive structure function KP, and can be written as the sum of 
several distinct contributions. In our fit we include the following four pieces: 

F? = Fgg + Ftfg + AF,£ + A*J. (13) 

Here the first and the second term, as indicated by the subscripts and by the superscripts, de- 
note the production of a quark-antiquark pair and the production of a quark-antiquark-gluon 
system with transversely-polarized photons. The third term takes into account the production of a 
quark-antiquark pair from a longitudinally-polarized photon, and the prefix A indicates that this 
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contribution belongs to higher twist (twist four). We have also included a transverse higher- twist 
contribution to qq production, denoted by AF^ 

Let us discuss these terms in more detail. To begin with the Q 2 dependence of (11), we reca- 
pitulate that in the leading-twist transverse contribution to qq production there is no log((5 2 /Qo)~ 
enhancement from the phase-space integral, whereas qqg production is of higher order in a s and has 
an a s ln(Q 2 /Qq) dependence. The third term, the longitudinal cross section of the qq final state, 
belongs to higher twist, and the phase-space integral provides a log(Q 2 /QQ) enhancement. The 
reason why this contribution is essential will be discussed below. The longitudinal contribution of 
the qqg final state is again of leading twist, but the logarithm is absent. It is therefore subleading 
in comparison with the transverse contribution, and will be disregarded in our parametrization. 

To get an estimate of the (3 spectrum, we consider the limits (3 — > 1 and j3 — > 0. Contact 
with the other variables is made through the kinematic relation (6). For the (3 dependence of 
the longitudinal cross section we can use (10) and (5): for the transverse case the situation is 
slightly more complicated, and we have to use both (9) and (11), in combination with (5). A 
more intuitive argument can be derived from the wave functions, and goes as follows. The limit 
M — > 0, which is the same as (3 — * 1, is related to the small-fcj behavior of the cross section. Let us 
return to the wave function (4). For both of the limits It — * and It — > oo, one finds that when kt 
approaches zero D^> vanishes with the same power in kt as the orginal wave function. (In eq.@, 
where the denominator is quadratic in kt, the integration over the azimuthal angle of It leads to 
the final cancellation of nonvanishing contributions.) This means that characteristic features of the 
light-cone wave functions (^ /0 , \E ,± , \E"^) remain unchanged after scattering. We expect that these 
results also hold for multi-gluon exchange. It is important to note that a single gluon (or photon) 
exchange, as opposed to the color-singlet two-gluon exchange, leads to a rather different spectrum 
at j3 ~ 1. This is because, instead of the second-order derivative in (4), in this case only the first 
derivative of the wave function is needed, which, unlike the second derivative, does not vanish when 
kt approaches zero. 

Applying these arguments to our wave functions we find, first for the the transverse quark- 
antiquark production cross section, that the cross section behaves like (1 — (3) ~ Qkt from 
(|l]) and M ~ kt from (6)), i.e., it vanishes when M becomes zero. For the second contribution 
with one gluon in the final state and with the tensor structure ~ k^kt — 2\kt\ 2 g^' 1 ', the cross 
section vanishes like (1 — (3) 2 . The subsequent integration over the quark- antiquark final state 
introduces a further suppression, leading to a (1 — /?) 3 behavior. Finally, for the wave function (|2|) 
of longitudinally-polarized photons, we find that *° ~ Q, i.e., the cross section goes to a constant 
different from zero: this means that near (3=1 the longitudinal cross section dominates and cannot 
be neglected. The other limit (3 — * or M — > oo is dictated by the high-energy behavior of the 
amplitudes, which is different for quark and gluon exchange. Spin-1/2 exchange is suppressed 
relative to spin-1 exchange, which leads to a dominance of gluon production at small (3 over the 
leading-order quark-antiquark production. We conclude that our three contributions, transverse 
qq and qqg production and longitudinal qq production, are important in rather distinct regions in 
(3, namely medium, small, and large (3, respectively. The transverse higher-twist contribution is 
expected to give a small negative correction which is due to phase space limitations at finite Q 2 pOfl . 

Finally, we comment again on the energy dependence. In contrast to the (3 spectrum, which 
can be traced back to rather general properties of the wave functions, perturbative QCD does not 
allow us to control the xp dependence of the cross section. In particular, n e ff for the leading- 
twist transverse cross sections cannot yet be predicted, and we therefore let the data decide on the 
preferred values of this exponent. As discussed at the end of the previous subsection, we expect 
a weak Q 2 dependence. For the higher-twist longitudinal part, on the other hand, theoretical 



9 



arguments have been given which indicate that the xjp dependence is given by the square of the 
gluon structure function at the momentum scale Q 2 /4/3, i.e., it should grow with Q 2 . In our fit we 
assume a universal xjp dependence of all higher-twist terms. The exponent is allowed to vary with 
Q 2 - 

After these remarks, we finally write down our Ansatz for the diffractive structure function. 
For the four terms in (11) we put: 



qg 



119 
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A 
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D 



Xjp 

\xjp 

XQ 

xjp 

XQ 

xjp 



II 2 
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112 
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Ql 
Q 2 

Q 2 



Ql 



In 



Q' 2 



ill 

Q 2 



^ In 



' Q 2 

MP 



+ 1.75 J /? 3 (1 - 2(3f 
+ 1.75 ) (3 3 {l-/3). 



(14) 



(15) 



The exponents are chosen to have the form 

m = n 2 o + n 2 i In 

ri4 = n4 o + ^4 i In 



In 


SH 


In 


SH 



(16) 



and a s is set to 0.25. 

As we have discussed before, the shapes of the /? spectra are restricted by properties of the 
light-cone wave functions, and one finds for the parameter 7 in the value 3. But, since the HI 
analysis in M reports a rather hard gluon distribution inside the Pomeron, we allow the parameter 7 
to deviate from the model prediction. For the leading-twist (transverse) cross sections the exponent 
?i2 will be left to the fit. We have introduced a simple Q 2 -dependent function, although a weak 
dependence is expected theoretically. The xjp dependence of the-higher twist contributions (724), on 
the other hand, is expected to be given by the square of the gluon structure function, and in principle 
we could use the measured gluon structure function. We have inserted an extra log(Q 2 / (4PQq)) in 
order to simulate the effect of having a structure function. For 714 we assume the same functional 
form as for but expect to find a stronger rise with Q 2 than for the leading-twist contributions. 
The parameter xq was introduced to minimize any effect from 712 and 714 on the overall Q 2 shape. 
The scale parameter Qq is taken to be lGeV 2 . 

Finally, we mention that, in principle, one should also allow for higher-twist corrections to qqg 
production. However, even if we assume a universal xjp dependence for all higher-twist pieces, this 
would increase the number of free parameters. We have found that, with the presently available 
amount and accuracy of data, it is not yet possible to determine these parameters with sufficient 
accuracy. Therefore, these higher-twist pieces will not be included in our fit. 



3 Fits to the Data 

ZEUS Data: We begin with a fit to the ZEUS data (l[, whose results are shown in Table 1 and 
Fig. 4. Some of the parameters have been restricted to a physically meaningful range by imposing 
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Figure 4: Fit compared to 1994 ZEUS data. Upper solid line: total result, dashed line: F^, dotted 



line: F^ g and dashed-dotted line: F^. 
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0.79 
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0.43 


0.12 





1 


xVd.o.f. 


12/43 



Table 1 



upper and lower limits. We note that the form of the qqg wave function in (3) leads to 7 = 3: our 
fit, cf. the value of 7 in Table 1 and the dotted curves in Fig. 4, indicates that the data prefer a 
larger value. The behavior of F q q g near (3 = 1 is therefore far from a 'hard gluon distribution' in 
the Pomeron. Our fit shows that it is possible to describe the ZEUS data without introducing a 
hard gluon inside the Pomeron. 
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Figure 5: The (3 spectrum at 
fixed xp = 0.001. Upper 
solid line: total result, dashed 
line: F^, dotted line: 
dashed-dotted line: 



F T - 
Fqq and 
lower solid line: AF^. 
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Fig.5 shows the (3 spectrum. One recognizes the subdivision into three distinct regions: the 
small-/? region with qqg production, the medium-/? region with transversely-produced qq pairs, and 
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the large-/? region where the longitudinal part dominates. The sum of all three contributions leads 
to a rather flat (3 spectrum. 




The Q 2 distribution in Fig. 6 has the expected shape: it is logarithmically increasing at j3 = 0.1 
(which we associate with qqg production), constant for the leading- twist transverse part at [3 = 
0.5, and decreasing for the longitudinal part at /3 = 0.9. The negative slope of the longitudinal 
contribution, which is due to its higher-twist nature, is partly compensated by the counteracting 
logarithm in Q 2 . The higher- twist correction for the transverse part, which is also decreasing when 
the absolute value is considered, leads to a small positive slope at (3 = 0.5 in the combined result. 

The transverse higher-twist correction to qq production wants to be negative, as expected from 
theoretical arguments [2C]. Its absolute magnitude is rather small, as can be seen in Figs. 5 and 
6, and it represents a correction to the three leading pieces. As mentioned before, we have also 
attempted to include higher-twist corrections to the qq production: with the available statistics it 
is not possible to assert whether such a contribution is present or not. 




Figure 7: Pomeron intercept 
for leading twist (solid line) 
and higher twist (dashed line) . 
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Although some of the parameters ri2o, \ have substantial errors, their values, as deter- 

mined from the fit, lead to an appealing theoretical scenario. The value ri2 o = 1-11 corresponds 
to an effective Pomeron intercept ajp(O) = 1.055, which is consistent with the 'soft' Pomeron of 
Donnachie and Landshoff. It exhibits, however, a slight rise with Q 2 , which suggests that in the 
leading-twist transverse cross section the Pomeron is already a mixture of the soft Pomeron and 
the gluon structure function. According to our discussion above, this means that the effective mo- 
mentum scale of the aligned-jet configuration in qq production is slightly higher than the typical 
hadronic scale of the soft Pomeron. The intercept for the higher-twist contribution (the dashed 
line in Fig. 7), on the other hand, shows a strong rise with increasing Q 2 . It tends towards a value 
of 1.2 at large Q 2 , which is compatible with the gluon structure function, i.e., the Pomeron is hard 
in this case. 

It is also interesting to compare the ZEUS fit with data from HI. Fig.8 shows all HI data points 
with the ZEUS fit result overlaid. The agreement is quite good when all points above xp = 0.01 
are ignored. The strongest deviation is observed in the (5 = 0.2 bin at low Q 2 . Here the ZEUS 
fit lies considerably above the HI data. The overall impression, however, is that in the kinematic 
range where the Pomeron dominates over secondary trajectories, both data sets are consistent. We 
find remarkable the good agreement between the HI data and the ZEUS fit at (3 = 0.9 and low Q 2 , 
since this region is not covered by ZEUS data. The fit seems to provide here a good extrapolation. 
In summary, we have shown that it is possible to describe the ZEUS data using the first three terms 
in (12),the fourth term providing a small correction. In particular, there is no need for a hard gluon 
inside the Pomeron. Furthermore, our fit is quite consistent with our theoretical expectations: the 
parameter 7 as well as the exponents n-2 and n.4 have chosen values which seem to confirm the ideas 
outlined in the previous section. 

HI Data: Next we turn to the HI data and describe the results of our fit to them. Most 
remarkable is the fact that we find two different solutions. The parameters of the first solution 
are given in Table 2, and the corresponding comparison of fit and data can be found in Fig.9. 
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139/130 



Table 2 



In order to avoid secondary trajectories, we have imposed an upper cut on xp of 0.01. A striking 
feature of this solution is the fact that the exponent 7 is well below 1, which can be interpreted 
as implying an initial gluon distribution that is singular for (3 = 1. The explanation for this result 
can be inferred from Fig.10, where we have compiled all relevant diagrams, i.e., the (3 spectrum, 
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Figure 8: The ZEUS fit compared to HI data. Upper solid line: total result, dashed line: F^, 
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Figure 9: The HI fit compared to HI data. Upper solid line: total result, dashed line: F^, dotted 



line: F^ g and dashed-dotted line: F^. 
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the Q 2 dependence and ap, in one figure. 





The HI data prefer a positive slope in Q 2 , even at a large 0, of 0.5. With the present parametriza- 
tion, this can only be achieved by making the gluon contribution (F^ g ) large. Of course, our rather 
simple approach is not competitive with the more sophisticated analysis that has been performed 
by HI Q. Nevertheless, it mimics the effect of an evolving singular gluon distribution rather well. 
The P spectrum in Fig. 10 shows the dominance of the qqg contribution over the transverse qq 
contribution, which also spreads into the large-/3 region. The longitudinal contribution has not 
disappeared completely, but is roughly a factor 2 smaller than in the fit to the ZEUS data. Since 
the qqg contribution has only a rather low Pomeron intercept associated with leading twist, the 
intercept for the higher-twist contribution is forced to a very high value of 1.3, in order to accom- 
modate the data. It is completely flat, because we demanded a positive slope in our fit. Without 
a lower limit on 7141, the slope would become negative. 

Comparing with the ZEUS data, we find again reasonable agreement between the two data sets. 
The most significant deviation is found for the lower Q 2 bins. We point out that, when one fits the 
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Figure 11: The HI fit compared to ZEUS data. Upper solid line: total result, dashed line: F^, 
dotted line: F^ g and dashed-dotted line: F^. 
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ZEUS data starting from the HI data as input, the parameter values always move back towards 
those found in the earlier fit, i.e., those in Table 1. In particular, the ZEUS data do not favour a 
value of 7 smaller than 1. If one requires a small value of 7, one obtains an acceptable fit only with 
an unreasonably small value for xq. If, in addition, xo is fixed at, say, 10~ 3 , the \ 2 value of the fit 
is substantially increased. 

As we have already said, the HI fit has a second local minimum for a large 7 of 8.5. The x 2 
is not much worse than in the previous fit (150 compared to 139). The parameters if this solution 
are given in Table 3, and the comparison with data is shown in Figs. 12 and 13. The (3 and Q 2 
distributions are shown in Fig. 14. On the whole, the distributions look similar to the ZEUS fit. In 
more detail, however, there are a few differences: the Q 2 shape for {3 = 0.5 is completely flat, and 
the Q 2 dependence of 77,4 is also flat. 
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x7d.o.f. 
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Table 3 



Summarizing the HI fit, we find when fitting our model to the HI data a solution which 
allows for the singular gluon interpretation, but there is also a second solution that is close to 
our model, i.e., consistent with the interpretation described in the context of the ZEUS data. 
Since the \ 2 values of both solutions are not that different from each other, one has to search for 
further consistency checks: more decisive tests might be provided by comparisons with the vector- 
meson production cross section or with hard diffractive jets. A singular gluon would give a strong 
transverse component, whereas in our model the longitudinal component is more pronounced. 
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Figure 12: The second HI fit, with the same notations as before. 
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Figure 13: The second HI fit compared to the ZEUS data. 
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4 Conclusions 

In this analysis of diffraction in deep-inelastic scattering at HERA, we have proposed a simple 
parametrization. Working in the wave-function formalism and starting from the 'hard part' of the 
diffractive cross section, we have suggested a simple extrapolation into the 'soft' nonperturbative 
region. The main feature of the parameterization is the decomposition of the (3 spectrum into three 
contributions which reside in separate regions with only little overlap: qqg production at low (3, 
transverse qq production at medium (3, and longitudinal qq production at large (3. These can be 
derived from the corresponding wave functions for the qqg and qq Fock states. The longitudinal 
contribution is higher twist, whereas the other two contributions are leading twist. The very 
different f3 dependence of these three contributions allows to determine their relative importance 
through the fits to the data. We fitted therefore our model to the recent ZEUS and HI data. 

An important result of the fit to the ZEUS data is the dominance of the Pomeron-quark coupling 
at larger (3 values {(3 > 0.3). The Pomeron-gluon coupling becomes substantial at lower (3 {(3 < 0.3). 
The region of large (3 ((3 > 0.9) is dominated by the longitudinal contribution. 
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The fit to the HI data leads to two solutions. In order to be open towards the HI conjecture of a 
singular gluon distribution at (3 = 1 H, we have allowed, in the second term of our parametrization, 
the exponent 7 to be variable. A singular gluon distribution would predict 7 < 1, whereas our 
model suggests 7 = 3, or even a little bit higher when Q 2 evolution is included. Whereas the ZEUS 
fit found only the option of a large 7 value, the two HI solutions have very different 7 values. 
In the first fit, 7 takes a small value, consistent with the HI conjecture. The Pomeron couples 
predominantly to gluons and only weakly to quarks. In the second solution, which has almost the 
same probability (x 2 value) as the first one, 7 is much larger and the Pomeron couples mostly to 
quarks. The coupling to gluons becomes substantial only at low f3, quite similar to the ZEUS fit. 

So far one cannot draw a final conclusion about the correct interpretation of both data sets. 
On the one hand, the conjecture of a singular gluon seems unlikely, in view of the ZEUS data. On 
the other hand, our parametrization with the two-gluon exchange as a model for the Pomeron can 
describe both data sets. More insight into the final state is needed, which could, for example, be 
provided by careful analyses of vector meson production and/or diffractively-produced jets. 
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